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HO 要 :结构 体 入 水 作为 一 个 高 度 非 线 性 的 冲击 动力 学 问题 ,广泛 存在 于 船 船 海洋 和 军事 等 工程 领 
域 ,因此 ,对 该 问题 进行 精确 而 高 效 的 数值 模拟 是 十 分 必要 。 采 用 有 限 质 点 法 (finite particle meth- 
od ,FPM) ,结合 无 反射 边界 ,数值 模拟 研究 了 二 维 结构 体 入 水 的 动态 过 程 。 通 过 分 析 模 形体 入 水 所 
受 的 垂 向 力 、 速度 和 压力 ,并 与 相关 文献 结果 对 比分 析 , 验 证 了 该 方法 的 有 效 性 。 在 此 基础 上 ,对 不 
同 工 况 下 的 入 水 问题 进行 了 模拟 ,研究 了 不 同 尺寸 .形状 、 初 速度 和 质量 等 情况 下 ,结构 入 水 冲击 时 
垂 向 力 和 速度 随时 间 的 变化 规律 ,以 及 不 同时 刻下 压力 响应 的 变化 趋势 。 研 究 结果 表明 : 枢 形 体 的 
底 升 角 尺寸 、 结 构 体 形状 、 入 水 初速 度 和 质量 等 因素 都 会 直接 影响 结构 体 入 水 的 力学 特性 ,为 入 水 
结构 体 的 力学 设计 提供 了 参考 。 
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Numerical simulation of water-entry based 
on the finite particle method 
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Abstract : Water entry of structures,a highly nonlinear impact dynamic process ,1s widely used in maritime 
and military fields. It 1s very necessary to simulate this problem accurately and effectively. The finite parti- 
cle method ( FPM) combined with the non-reflection boundary was adopted to simulate the problem of wa- 
ter entry. The vertical force, the velocity and the pressure of two-dimensional wedge were analyzed and 
compared with the results of the literature to verify the effectivity of this method. The factors for the prob- 
lem of water entry were also studied in this paper. The effects of size, shape , initial velocity and mass on 
vertical force , velocity and pressure at different time were analyzed. It is shown that the mechanical proper- 
ties of the structure are affected by the bottom angle of the wedge , the shape of the structure , the initial ve- 
locity and the mass , which has a reference for the design of the water-entry structure. 
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入 水 问题 作为 一 个 常见 的 物理 现象 ,广泛 地 存 
在 于 工程 工业 和 生产 生活 中 ,尤其 是 在 水 上 交通 和 
军事 等 方面 ,如 轮船 ,水 上 飞机 、 潜 舰 导 弹 入 水 等 问题 ， 
都 会 涉及 到 结构 体 运动 与 水 面 的 撞击 ,由 于 结构 体 受 
到 较 大 的 水 动力 作用 ,会 发 生 失 稳 、 变 形态 至 破损 的 现 
象 ,从 而 直接 影响 入 水 设施 及 系统 的 安全 稳定 运行 。 

结构 体 和 水 是 一 个 复杂 的 、 典 型 的 流 固 炮 合 问 
题 ,在 入 水 过 程 中 ,结构 体会 对 水 体 产 生 冲 击 作用 ， 
造成 水 的 飞溅 和 流动 ,同时 水 体 又 会 产生 反作用 力 
施加 在 结构 体 上 ,使 之 发 生变 形 或 者 破坏 … 。 自 文 
献 [2] 和 文献 [3 对 入 水 问题 展开 人 研究 以 来 ,国内 外 
学 者 采用 不 同 的 方法 对 该 问题 进行 了 深入 分 析 “。 
文献 [4] 对 横 形 体 跌 落 入 水 过 程 进行 实验 人 研究 ,分 
析 了 结构 体 上 的 撞击 载 集 ,并 与 二 维 数值 模拟 结果 
进 每 了 对 比 。 文 献 [6] 通 过 实验 和 数值 模拟 的 方法 
人 研 放 了 弹性 棉 形 体 的 入 水 冲击 问题 。 文 献 [7] 对 细 
长 体高 速 入 水 物体 的 表面 压力 特性 进行 了 试验 研究 。 

随 关 数值 计算 方法 的 快速 发 展 ,数值 模拟 也 成 
为 研究 人 水 问题 的 重要 手段 之 一 ”。 由 于 结构 体 的 
he 形状 和 材料 等 参数 会 直接 影响 人 水 的 动态 过 
程 Y 研 究 人 员 采 用 不 同 的 模拟 方法 对 不 同 工 况 下 的 
人 水 问题 进行 了 研究 “”。 基 于 网 格 法 的 数值 模拟 方 
法 车 被 广泛 采用 ,文献 [10j] 应 用 动 网 格 技 术 和 流体 
体积 (volume of fluid, VOF) 方 法 模拟 了 平底 结构 日 
而 入 水 及 空 泡 演化 过 程 。 

近年 来 ,无 网 格 方法 也 被 应 用 于 入 水 问题 的 研 
务 s= 文 献 [ 11 ] 采 用 任意 拉 格 明日 - 欧 拉 算法 ALE 
(.Xrbitrary Lagrange-Euler) 模 拟 了 二 维 棉 形 体 入 水 过 
Tes" SCR 12-13 |] 则 采用 光滑 粒子 流体 动力 学 方法 
(SPH) 研究 了 该 问题 。 文 献 [14] 采 用 Godunov SPH 
方法 模拟 研究 了 弹性 体 的 入 水 过 程 。 针 对 传统 SPH 
方法 计算 精度 低 、 边 界 区 域 计 算 不 足 .数值 不 稳定 等 
问题 ,LIU 等 ”提出 了 一 种 精度 更 高 的 改进 方 
法 , 即 有 限 质 点 方法 (finite particle method ,FPM ) , 适 
合 于 模拟 具有 有 自由 面 移 动 界面 和 可 变形 边界 的 大 变 
形 流动 问题 “” 。 文献 [17] 采 用 FPM 方法 研究 了 孤 
立波 的 扑 坡 过 程 ;文献 [19-20 ] 分 别 采 用 FPM 方法 研 
究 了 水 中 结构 体 与 波浪 的 相互 作用 。 然 而 ,目前 采用 
FPM 方法 模拟 入 水 问题 的 相关 研究 报道 较 少 。 

本 研究 采用 FPM 方法 对 入 水 问题 进行 了 数值 
模拟 ,通过 结合 无 反射 边界 消除 了 水 箱 壁 面 回 弹 波 
的 影响 ,与 已 有 研究 结 采 进行 对 比 , 验 证 了 该 方法 的 
有 效 性 和 准确 性 ;深入 探究 了 结构 体 的 底 升 角 、 形 


状 、 初 速度 和 质量 等 因素 对 入 水 过 程 的 影响 ,为 入 水 
冲击 等 复杂 工程 问题 的 设计 优化 提供 理论 依据 。 


1 FPM 方法 


1.1 FPM 基本 理论 


根据 文献 [16] 提出 的 FPM 方法 的 基本 理论 ,二 
维 坐标 系 下 FPM 基本 方程 的 离散 格式 为 5221 
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(3) 
式 中 :大 (nm) 为 粒子 i 的 函数 及 其 导数 ;f(r) 为 粒子 
j FS PR ZR Wye DN RR C i SG. f(xi,y;) ， 
fi Gy) 和 (zwi,Yi) 分 别 为 点 (x, y). 处 的 函数 及 
x Ay Jr 18] E FS fi S306; m; 和 pj 分 别 为 粒子 j 的 质量 
和 密度 ; N 为 文 持 域内 与 粒子 相互 作 用 的 粒子 7 的 
总 数 ;4 为 修正 和 矩阵; a I B 代表 坐标 方向 ,以 a 为 


例 , 当 waw = 0 时 表示 困 数 本 有 身 , 当 aw = x 时 ,分 别 表 
ZI RAE x y 方向 上 的 导数 。 
本 研究 中 采用 了 应 用 广泛 的 B 样 条 函数 
2 1 
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式 中 :二 维 空间 中 e, 为 13MTmj ; h 为 光滑 长 度 ; 
R = r/h ,7 为 两 个 粒子 之 间 的 距离 。 
1.2 ”控制 方程 的 FPM 离散 形式 

对 于 流 固 耦 合 问题 而 言 , 由 于 粒子 的 不 一 致 性 ， 
传统 SPH 方法 存在 精度 不 高 的 问题 ,文献 [16 ] 


担 岂 的 FPM 方法 ,用 于 改进 SPH 方法 固有 的 粒子 不 
- 右 性 ,提高 了 求解 精度 。Navier-Stokes( N-S) 方程 


的 SRPM 离散 形式 为 
e m 
di — -P ox" (5) 
du __ p , ðe” 
P di Ox" d ax! 


CD 为 了 克服 压力 的 不 稳定 性 ,在 连续 性 方程 和 动 
量 年 恒 方 程 中 分 别 加 入 密度 扩散 项 和 人 工 儿 性 项 ， 
RE) 可 以 改写 为 
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xu 为 速度 矢量 ,wy =u; -u ,hy = (s € h)72, 
y; =2(p; -p))xi/ |r; 1” ;扩散 项 系数 6 20.1; ATA 


zi 为 粒子 的 空间 坐标 位 置 ; po 为 初始 密度 ; c 为 
声速 。 

本 研究 将 流体 视 为 弱 可 压缩 流体 ,其 压力 为 “ 

p -c(p-p) (7) 


1.3 无 反射 边界 处 理 


结构 体 入 水 冲击 时 会 产生 冲击 波 , 冲击波 到 达 
箱 体 边界 后 会 回 弹 并 干扰 内 部 流 场 ,从 而 对 结构 体 
的 运动 产生 一 定 的 影响 ”。 为 了 尽 可 能 避免 上 述 
干扰 ,采用 粒子 法 模拟 时 通常 需要 水 体 的 尺寸 远大 
于 结构 体 的 尺寸 ,使 得 计算 区 域 较 大 ,导致 计算 成 本 
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增加 。 本 人 研究 在 水 体 边界 选取 一 定 厚度 的 粒子 层 ， 
处 理 为 无 反射 边界 ,从 而 可 以 完全 避免 水 箱 壁 面 边 
界 导 致 的 冲击 波 反 射 干扰 现象 ,能够 减 小 水 体 的 尺 
寸 并 大 幅度 提高 计算 效率 ”。 经 多 次 实验 模拟 , 结 
果 表 明 , 选 取 的 粒子 层 厚 度 为 搜索 区 域 的 10 倍 为 
IgE. 
在 所 选 的 边界 层 内 ,动量 守恒 方程 为 
du — 1 (ge) (8) 


AP: n =0.5, C = 05,K, 210x105, 


2 数值 计算 验证 


2.1 算 例 模型 的 介绍 


为 了 检验 本 人 研究 FPM 方法 的 准确 性 ,对 模 形 体 
和信 水 问题 进行 了 模拟 研究 , 算 例 模型 如 图 Ir. 
水 体 的 宽度 为 2 m, 深度 为 0.7 m, 初 始 密 度 为 
1 000 kg/m? ;将 析 形 体 处 理 为 刚性 结构 ,其 初速 度 为 
6. 15 m/s, JEJ a 为 30°, 质 量 为 241 kg。 粒 子 直 
径 d 为 0.002 m, 共 有 359 748 个 粒子 ,hh 为 1.23 d,c 
为 200。 为 保证 求解 的 稳定 性 ,时 间 步 长 应 满足 
CFL 条 件 , 即 At<0.5 x h/c, At BO 5 x10 s RER 
计算 的 总 时 长 为 0.025 s, 


(a) 入 水 模型 


B=500 mm 


(b) BUE HE BUR 


图 1 算 例 模型 的 示意 图 


Fig.1 The schematic diagram of the computation model 
2.2. 算 例 模 型 的 验证 


本 研究 采用 FPM 方法 模拟 研究 了 模 形 体 撞击 
入 水 的 动态 过 程 ,并 与 文献 [4 |] 和 文献 [12] 进 行 了 
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对 比分 析 。 图 2 为 作用 于 模 形 体 的 垂 向 力 曲 线 , 图 3 
给 出 了 不 同时 刻下 横 形 体 斜 边 的 压力 分 布 曲 线 。 图 
中 的 压力 分 布 按 上 述 文献 进行 了 相应 的 无 量 纲 处 
理 , 其 中 , F, 为 攀 形 体 所 受阻 力 和 重力 之 差 , 即 垂 向 
JJ. y 为 枢 形 体 斜 边 纵 坐 标 , y, 为 枫 形 体 顶 点 的 纵 坐 
标 , ya 为 棉 形 体 底 边 上 的 高 ,h 为 水 面 高 度 , S 为 水 面 
高 度 与 棉 形 体 项 点 纵 坐 标 之 差 ,V 为 结构 体 的 速度 。 


"= ZHAO 实验 值 


9 000[ 一 一 ZHAO 理论 值 

g 000[--- GONG SPH-W™* 结 果 (4=0.001 2 m) 
c— SPH-KGC 结 果 (4d=0.002 m) 

7 o0ol-—- FPM 结 果 (q=0.002 m) 


0 
0.000 


0.005 0.010 0.015 


t/s 
图 2 枫 形 体 的 垂 向 力 随 时 间 的 变化 曲线 


Fig.2 Variation of the vertical force on the wedge 


0.020 0.025 


HE2 可 知 , 枫 形 体 所 受 的 垂 向 力 随时 间 的 变化 
趋 攻 与 文献 结果 较为 吻合 。 根 据 垂 癌 力 的 变化 趋势 ， 
将 过 程 分 为 3 个 阶段 :入 水 初期 (0<t <0. 008 33 s) ,年 
向 汶 急 剧 增 大 ;入 水 中 期 (0.008 33 s: «0.016 67 s), 
重 岂 力 达 到 顶峰 并 出 现下 降 趋势 ; 入 水 后 期 
(07016 67 s1z0. 025 s) ,3f [8] 7] AHIR N 
为 了 进一步 验证 模拟 结果 ,本 人 研究 选取 3 个 阶 
段 串 的 不 同时 刻 , 给 出 了 沿 着 横 形 体 斜 边 的 压力 分 
dpi] 3 所 示 。 
图 2 和 图 3 还 给 出 了 理论 解 和 基于 SPH 方法 改 
进 的 SPH-W"” (moving-least-squares kernel ) 方法 和 
SPH-KGC ( kernel-gradient-correction ) 方法 模拟 所 得 
结 采 。 从 图 2 可 以 看 出 ,不 同方 法 所 得 的 垂 回力 与 
文献 14] 的 理论 解 及 实验 值 基本 吻合 ,但 在 时 间 步 
又 0.008 ~0. 017s 的 范围 内 均 大 于 实验 结 末 。 观 察 
不 同时 刻 的 压力 分 布 曲线 图 3 可 以 看 出 :在 入 水 初 
期 ,不 同方 法 得 到 的 压力 峰值 之 间 的 差异 较 大 ;在 和 人 
水 中 期 ,本 研究 FPM 方法 和 文献 [12] 的 SPH-W 
方法 与 理论 解 的 压力 分 布 趋势 基 本 一 致 ,这 2 种 方 
法 所 得 结果 差异 较 小 ,与 理论 解 较 为 吻合 ,而 FPM 
方法 所 得 曲线 波动 小 ,较为 光滑 平缓 。 此 外 ,FPM 方 
法 比 SPH-W ”方法 较 接 近 于 实验 值 ,虽然 SPH-KGC 
方法 最 接近 实验 值 ,但 与 理论 值 偏 离 很 大 。 在 入 水 
后 期 , SPH-W “方法 所 得 压力 曲线 波动 较 大 , 而 
FPM 方法 的 曲线 波动 小 ,与 理论 解 相 吻合 ,也 比 


SPH-W ”方法 接近 实验 值 。 同 样 SPH-KGC 方法 虽 
然 最 接近 实验 值 , 但 仍然 与 理论 值 偏离 很 大 。 
14 : 
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图 3 不 同时 刻下 横 形 体 斜 边 的 压力 分 布 曲线 
Fig.3 Pressure distribution on the wedge bevel 
本 研究 还 给 出 了 该 3 种 时 刻下 水 体 的 压力 分 布 
云图 ,如 图 4 所 示 。 从 图 4 可 以 看 出 ,采用 不 同 粒 子 
尺寸 得 到 的 压力 分 布 云 图 在 不 同时 刻下 的 变化 趋势 
几乎 一 致 , 横 形 体 对 水 体 的 冲击 作用 随 者 入 水 深度 
的 增加 而 减 小 ,这 是 因为 随 着 入 水 时 间 的 增加 ,在 水 
体 对 模 形 体 的 阻力 作用 下 , 模 形 体 的 动能 不 断 消 耗 
并 减 小 , 模 形 体 对 水 体 的 冲击 作用 逐渐 减 小 。 
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€ 图 4 不 同 粒子 尺寸 下 水 体 的 压力 分 布 云 图 
O Fig.4 Pressure distribution of water with different particle sizes 


为 了 分 析 粒 子 尺 寸 对 计算 结果 的 影响 ,这 里 还 
给 出 了 作用 于 枫 形 体 的 垂 向 力 曲线 和 不 同时 刻下 枫 
形体 斜 边 的 压力 分 布 曲线 ,分 别 如 图 5 和 图 6 所 示 。 

从 图 5 可 以 看 出 ,粒子 尺寸 对 槐 形体 的 垂 向 力 
影响 较 小 。 从 图 6 可 以 看 出 ,粒子 尺寸 对 横 形 体 的 
压力 峰值 及 其 出 现 的 时 间 有 一 定 影响 。 由 于 结构 体 
所 受 压力 主要 来 自 于 水 体 的 作用 , 随 着 粒子 尺寸 的 
增 大 ,构成 溅 起 水 花 的 粒子 个 数 越 少 ,尤其 是 入 水 初 
期 沿 着 斜 边 候 升 的 水 体 粒 子 较 少 ,从 而 出 现 支持 域 
内 粒子 数目 过 少 导 致 的 求解 偏差 。 由 图 6 可 知 , 当 
粒子 尺寸 等 于 0. 004 m 时 ,模拟 结果 与 理论 解 存在 
着 较 大 的 差异 。 随 着 粒子 尺寸 的 减 小 , 当 尺 寸 
为 0.002 m 时 ,除了 在 人 水 初期 会 出 现 压 力 峰 值 高 
于 理论 值 之 外 ,其 余 模 拟 结果 均 与 理论 解 较 为 吻合 。 
粒子 直径 为 0.001 5625 m 和 0. 002 m 在 i5-8250U 


的 CPU 上 计算 总 时 长 分 别 为 5.5h 和 3 h, 因 此 , 综 
合 考 虑 计算 效率 ,本 人 研究 选取 的 粒子 尺寸 为 0.002 m, 
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图 5 不 同 粒子 尺寸 下 横 形 体 的 垂 向 力 随时 间 的 变化 曲线 


Fig.5 Variation of the vertical force on the wedge 


with different particle sizes 
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图 6 不 同 粒子 尺寸 下 横 形 体 斜 边 的 压力 分 布 曲线 


Fig.6 Pressure distribution on the wedge bevel with different particle sizes 
行 了 分 析 。 
类 \ 4 
UK: 吉 构 体 的 尺寸 .形状 、 初 速度 和 材料 都 会 下 图 7 ~ 图 9 给 出 了 20°、25°、30°* 和 35° 4 种 不 同 


接 影响 入 水 的 动态 过 程 ,为 了 深入 人 研究 结构 体 和 水 ” 底 升 角 的 横 形 体 入 水 的 模拟 结果 ,包括 水 体 的 压力 
相互 作用 ,本 研究 模拟 了 不 同 工 况 下 不 同 结构 。” 分 布 云 图 、 垂 向 力 、 速 度 以 及 横 形 体 斜 边 的 压力 随时 
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体 的 入 水 过 程 ,并 对 结构 体 的 力学 特性 和 冲击 作用 ” 间 变 化 的 曲线 图 。 
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图 7 不 同 底 升 角 下 水 体 的 压力 分 布 云图 


Fig.7 Pressure distribution of water with different bottom angles 


其 他 参数 与 2.1 中 的 数值 保持 一 致 。 由 图 7 可 
知 , 模 形体 入 水 后 ,水 体 开 始 沿 斜 边 爬 升 , 模 形 体 两 
侧 压 力 分 布 呈现 出 较 好 的 对 称 性 , 随 着 人 水 时 间 的 
增加 , 模 形 体 神 击 水 体 产 生 的 压力 波 逐 渐 回 四 周 传 
TE BUE BUS Tff BU] BUE MOSEZK MSIE BIE 
S Ue le. , IKERRI B7 Hi HE BUD, Vols EE 
越 大 且 飞 溅 射流 越 细 长 。 

从 图 8 和 图 9 可 以 看 出 模 形 体 的 底 升 角 对 垂 向 
力 、 速 度 和 模 形 体 斜 边 所 受 压力 有 明显 的 影响 。 由 
图 8 可 知 , 棉 形体 的 垂 向 力 和 下 降 速 度 密 切 相 关 。 
随 着 横 形 体 底 升 角 的 减 小 , 横 形 体 所 受 的 垂 向 力 的 
峰值 越 大 且 越 早 达到 峰值 ,在 垂 向 力 的 作用 下 模 形 
体 的 下 降 速 度 随 着 底 升 角 的 减 小 而 减 小 ; 随 着 入 水 
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时 间 的 增加 ,水 体 对 槐 形体 的 垂 癌 力 逐 渐 减 小 后 又 
增加 , 模 形 体 速度 的 降幅 则 逐渐 增 大 而 后 减 小 。 

由 图 9 可 以 看 出 , 枢 形 体 底 升 角 越 小 , 随 着 人 水 
时 间 的 增加 , 棉 形 体 佟 边 上 的 压力 变化 越 大 。 在 入 
水 初期 , 析 形 体 斜 边 所 受 压力 及 其 峰值 随 着 抵 升 角 
的 减 小 而 增 大 ,但 在 入 水 的 中 期 和 后 期 ,的 升 角 为 
20° 的 横 形 体 和 斜 边 所 受 压力 最 小 。 由 图 9 可 知 ,在 和 人 
水 冲击 的 中 斯 和 后 期 , 棉 形 体 斜 边 上 的 压力 呈现 出 
高 频 的 震荡 ,这 是 因为 棉 形 体 冲击 过 程 中 形成 压力 
波 造 成 水 体 发 生 波浪 起 伏 ,从 而 使 得 棉 形 体 斜 边 所 
受 压 力 出 现 波 动 。 由 此 可 以 说 明 , 模 形体 入 水 冲击 
所 受 垂 向 力 、 速 度 和 压力 部 与 其 抬升 角 具有 紧密 的 
关联 。 
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图 8 不 同 底 升 角 下 槐 形体 的 垂 向 力 和 速度 随时 间 的 变化 曲线 


Fig.8 Variation of the vertical force and velocity on the wedge with different bottom angles 
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图 9 不 同 底 升 角 下 棉 形 体 斜 边 的 压力 分 布 曲线 


Fig.9 Pressure distribution on the wedge bevel with different bottom angles 


影响 

由 于 入 水 产生 的 水 溅 、 垂 回力 和 速度 都 会 随 着 
结构 体形 状 的 改变 而 产生 较 大 差异 ,因此 ,本 人 研究 选 
取 了 圆柱 体 、 柜 形体 和 立方 体 3 种 入 水 结构 ,其 横 截 


3.2 形状 


面 的 尺寸 分 别 为 : 圆 形 直径 为 0.2 m, 三 角形 斜 边 长 
为 0.2m, 底 升 角 为 30? ,方形 边 长 为 0.2 m。 结 构 体 
的 质量 均 为 50 kg, 和 人 入水 初速 度 为 6.5 m/s, 

从 图 10 的 压力 分 布 云 图 可 以 看 出 ,立方 体 入 水 
过 程 中 与 水 体 的 冲击 接触 面 较 大 ,导致 水 体 受到 了 
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较 大 的 冲击 作用 ,冲击 波 向 抵 部 水 层 传 播 的 同时 问 
两 边 移 动 ,因此 水 体 中 出 现 较 为 明显 的 震荡 弧 形 压 
力 波 。 横 形体 和 圆柱 体 入 水 过 程 中 水 体 的 压力 主要 


P: 0 12000 24 000 36 000 48 000 60 000 72 000 84 000 P: 


P: 0 12000 24 000 36 000 48 000 60 000 72 000 84 000 P: 


0 12000 24 000 36 000 48 000 60 000 72 000 84 000 P: 


0 10000 20 000 30 000 40 000 50 000 60 000 70 000 0 4000 8000 12000 16 000 20 000 24 000 28 000 


Ce Cl 


0 10000 20 000 30 000 40 000 50 000 60 000 70 000 P: 0 4000 8000 12000 16 000 20 000 24 000 28 000 


ChinaXiv 合 作 期 刊 


集中 在 结构 体 底部 ,圆柱 体 产生 的 冲击 压力 略 小 于 
BUE MS ,造成 的 水 溅 也 较 小 ,而 模 形 体 由 于 其 目 身 的 
结构 特点 ， 
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水 体 沿 着 斜 边 更 容易 向 两 边 飞 溅 。 
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mu 图 10 不 同 结构 体 的 水 体 压 力 分 布 图 

is Fig. 10 The pressure distribution of water with different structures 


泛 图 11 和 图 12 分 别 给 出 了 不 同 结构 体征 向 力 、 
a RIJG 图 中 ,x 为 底部 中 心 
点 湛 边 缘 的 结构 半 宽 。 从 图 11(a) 可 知 :立方 体 的 
证 河 力 整体 较 小 ;圆柱 体 在 入 水 瞬间 垂 向 力 最 大 , 随 
后 不 断 震 荡 减 小 ;机 形体 的 垂 向 力 则 先 增 大 后 缓慢 
减 小 。 结 合 图 11(b) 可 以 看 出 :由 于 立方 体 所 受 的 
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向 上 的 垂 疝 力 最 小 ,其 下 降 速 度 不 断 增加 ;圆柱 体 和 
棉 形 体 的 速度 逐渐 减 小 ,其 中 ,后 者 入 水 速度 的 变化 
幅度 较 大 ,这 是 因为 在 冲击 作用 下 水 体 沿 着 其 斜 边 
问 两 侧 流动 ,二 者 的 接触 作用 面积 增 大 ,加 大 了 水 体 
对 横 形 体 的 阻力 ,导致 其 向 下 的 速度 最 小 。 
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图 11 不 同 结构 体 的 垂 向 力 和 速度 随时 间 的 变化 曲线 


Fig. 11 


Variation of the vertical force and velocity with different structures 


由 图 12 可 知 ,由 于 水 体 的 反 向 冲击 作用 ,结构 
体 在 与 水 体 接触 的 区 域 具有 较 大 的 压力 。 在 入 水 过 
程 中 ,圆柱 体 所 受 的 压力 整体 变化 较 小 ,中 心底 部 的 
压力 最 大 , 沿 弧 两 边 逐 渐 减 小 ;立方 体 所 受 压 力 均 匀 
分 布 在 与 水 体 接触 的 底面 上 ,并 且 整 体 随 着 入 水 时 
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间 的 增加 略微 增 大 。 对 于 横 形 体 而 言 ,入 水 初期 压 
力 最 大 , 约 为 其 余 2 种 结构 体 压力 的 2 倍 ,并 且 压 力 
沿 着 其 斜 边 逐 渐 增 大 ;随后 ,压力 大 幅 下 降 , 当 和 斜 边 
已 完全 入 水 后 , 冬 边 压力 分 布 较为 均匀 。 
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Fig. 12 Pressure distribution with different structures 


3:3- 初速 度 的 影响 


sq 为 了 分 析 入 水 初速 度 对 冲击 过 程 的 影响 ,本 研 
完 选 取 直 径 为 0.2m 的 圆柱 体 作 为 和 水 结构 体 ,对 
看 同 初速 度 下 的 入 水 过 程 进 行 了 模拟 研究 。 结 构 体 
的 帝 量 为 50 kg ,初速 度 Vo 介 于 2.5 ~7.5 m/s 之 间 。 
结构 体 的 垂 向 力 、 速 度 和 压力 随时 间 的 变化 曲线 如 
LIS 和 图 14 所 示 。 为 了 更 好 地 分 析 速 度 的 变化 规 
律 ,图 13(b) 中 给 出 了 结构 体 的 速度 相对 变化 量 。 

一 观察 图 13 可 以 得 出 ,初速 度 对 结构 体 的 垂 向 力 
和 速度 的 变化 趋势 影响 不 大 ,但 对 二 者 的 数值 有 一 
定 的 影响 。 随 着 初速 度 的 减 小 , 垂 回 力 越 小 ,速度 相 
对 变化 量 曲 线 越 平缓 ,这 说 明 初 速度 越 小 ,在 同一 时 
刻下 结构 体 的 冲击 越 小 ,水 体 对 结构 体 的 垂 向 力 越 
小 ,结构 体 的 速度 变化 越 小 ;由 于 初速 度 越 大 惯性 越 
大 ,当初 速度 大 于 等 于 5.5 m/s 时 , 垂 向 力 会 在 := 
0.0002s 附近 产生 一 个 峰值 且 啊 应 时 间 较 短 ,峰值 
随 春 初 速度 的 增 大 而 增 大 。 因 此 , 右 要 减 小 结构 体 
的 垂 回力 ,可 以 降低 结构 体 的 入 水 初速 度 。 

从 图 14 所 示 的 结构 体 的 压力 曲线 图 可 以 看 出 ， 
入 水 速度 是 影响 结构 体 承受 压力 的 主要 因素 之 一 。 
结构 体 的 入 水 初速 度 越 小 ,所 承受 的 压力 与 冲击 动 
能 的 比值 越 大 ; 束 压 力 的 整体 分 布 而 言 ,圆柱 体 入 水 
底部 的 压力 最 大 , 沿 着 弧 边 压力 逐渐 减 小 ,但 在 水 花 
溅 起 的 接触 区 域内 ,结构 体 表 面 压力 出 现 了 剧烈 的 
波动 ,甚至 出 现 了 较 大 的 压力 值 ,初速 度 越 小 ,波动 
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Fig.13 Variation of the vertical force and velocity on the 
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图 14 不 同 初速 度 下 圆柱 体 的 压力 分 布 曲线 
Fig. 14 Pressure distribution on the cylinder 


with different initial velocities 
3.4 质量 的 影响 


为 了 人 研究 结构 体质 量 对 冲 击 过 程 的 影响 ,本 研 
究 也 对 不 同 质量 下 的 结构 体 和 人 水 过 程 进行 了 模拟 研 
究 。 选 取 直 径 为 0.2m, 初 速度 为 6.5 m/s 的 圆柱 体 
作为 人 水 结构 体 ,结构 体 的 质量 介 于 10 ~ 50 kg 
Zi, 


~ d As。 / fr TD] ' 
RY nr NAIL- HI Gal 
NNA IVa | FRHT! 


图 15 和 图 16 分 别 给 出 了 结构 体 的 垂 癌 力 、 速 
度 和 压力 随时 间 的 变化 曲线 。 
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图 15 不 同 质量 下 圆柱 体 的 垂 向 力 
和 速度 随时 间 的 变化 曲线 
Fig.15 Variation of the vertical force and velocity 
on the cylinder with different masses 
由 图 15(a) 可 知 , 质 量 对 结构 体 的 垂 向 力 变化 
趋势 影响 不 大 ,结构 体 垂 向 力 的 最 大 值 和 二 次 峰值 
以 及 其 出 现 的 时 间 均 保持 一 致 ,在 人 水 中 期 和 后 期 
垂 回 力 随 质 量 的 增 大 而 增 大 。 由 图 15 (b) 可 以 看 
出 , 随 着 质量 的 增加 ,速度 曲线 的 变化 越 平 绥 , 这 说 
明 质 量 越 大 惯性 越 大 ,因此 冲击 速度 的 变化 较 小 。 
图 16 的 压力 特性 曲线 显示 ,质量 越 小 , 入水 初期 结 
构 体 底部 的 压力 越 小 , 随 着 和 人 水 座 度 增加 ,水 体 震 荡 
并 溅 起 水 花 , 质 量 越 小 的 结构 体 则 越 容易 受到 水 体 
震 沪 的 影响 ,从 而 导致 压力 出 现 剧烈 的 波动 ,而 质量 
较 大 的 结构 体 的 压力 曲线 则 相对 平缓 。 与 人 水 初速 
度 相 比 ,质量 对 结构 体 所 承受 冲击 压力 的 影响 相对 
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图 16 不 同 质量 下 圆柱 体 的 压力 分 布 曲线 


Fig. 16 Pressure distribution on the cylinder 


with different masses 
4 结 it 


本 人 研究 采用 FPM 方法 结合 无 反射 边界 处 理 方 
法 ,对 二 维 结构 体 入 水 问题 进行 了 模拟 研究 ,着 重 分 
析 了 结构 体 的 尺寸 .形状 、 入 水 初速 度 和 质量 等 因素 
对 入 水 冲击 过 程 的 影响 。 通 过 分 析 结 构 体 的 受 力 情 
况 和 速度 的 动态 变化 过 程 ,总 结 了 不 同 工 况 下 和 人 水 
过 程 的 运动 变化 规律 ,对 入 水 结构 体 的 外 形 设计 和 
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人 和信 水 参数 选取 提供 技术 支撑 ,对 理解 结构 体 和 水 体 
的 相互 作用 以 及 结构 体 的 优化 设计 提供 重要 参考 价 
值 ,所 得 结论 如 下 。 

1 ) 通 过 模拟 经 典 的 枫 形 体 和 人 水 过 程 , 与 相关 理 
论 解 和 现 有 SPH-W"^ .SPH-KGC 方法 进行 比较 ,发 
现 采 用 FPM 方法 结合 无 反射 边界 得 到 的 棉 形 体 斜 
边 压 力 和 垂 向 力 与 理论 值 具有 较 好 的 吻合 性 ,同时 
比较 接近 实验 值 。 从 理论 和 实验 层面 验证 了 本 研究 
FPM 方法 和 所 建 模 型 的 准确 性 和 有 效 性 。 

2 ) 对 横 形 体 的 底 升 角 进 行 了 分 析 , 随 着 底 升 角 
的 减 小 , 棉 形 体 佟 边 上 的 压力 变化 越 大 , 底 升 角 越 
小 ,压力 峰值 越 大 且 峰 值 出 现 的 越 早 ; 当 人 水 后 期 枫 
形体 稳定 后 , 底 升 角 越 小 ,结构 体 底面 的 压力 越 小 ; 
模 形 体 底 升 角 的 大 小 与 垂 向 力 和 速度 下 降幅 度 呈 和 负 
相关 。 

3) 结构 体 的 形状 对 入 水 结构 体 的 垂 向 力 、 速 度 
和 压力 均 有 较 大 影响 。 随 着 和 人 水 冲击 时 间 的 增加 ， 
圆柱 体 的 垂 向 力 震 荡 减 小 , 模 形 体 的 垂 回 力 先 增 大 
后 缓慢 减 小 ,而 立方 体 的 垂 向 力 整体 较 小 ; 圆柱 体 和 
横 形 体 的 入 水 速度 不 断 减 小 ,而 立方 体 的 速度 逐渐 
增 大 。 由 于 横 形 体 更 容易 造成 水 体 飞溅 ,在 入水 过 
程 中 结构 体 所 受 压 力 明 显 减 小 ;圆柱 体 对 水 体 的 挤 
压 作 用 较为 明显 ,因此 底部 所 受 的 冲击 压力 较 大 且 
降幅 较 小 ,压力 沿 弧 形 向 两 侧 区 域 不 断 减 小 ;立方 体 
的 平底 结构 对 压力 的 变化 趋势 无 明显 影响 。 

4) 结 构 体 的 初速 度 越 大 , 垂 回力 和 速度 相对 变 
化 量 的 降幅 均 增 大 ;结构 体 所 受 压 力 随 着 初速 度 的 
增加 相对 减 小 ,在 一 定 的 范围 内 , 增 大 初速 度 可 以 降 
低压 力 的 波动 震荡 ;结构 体 的 质量 越 大 ,水 体 对 结构 
体 的 减速 作用 越 小 , 垂 向 力 .速度 和 表面 压力 峰值 整 
体 均 随 大 质量 的 增 大 而 增 大 ;与 初速 度 相 比 ,结构 体 
质量 对 入 水 过 程 的 影响 相对 较 小 。 
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